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ABSTRACT 



The effects of grazing on soil organic carbon (SOC) dynamics, particularly in the 
tropics, are still poorly understood. Plant compensation to grazing, whereby plants 
maintain leaf area (C input capacity) despite consumption (C removal) by grazers, 
has been demonstrated in tropical grasslands but its influence on SOC is largely un- 
explored. Here, the effect of grazing on plant leaf area index (LAI) was measured in 
a field experiment in Serengeti National Park, Tanzania. LAI changed little for graz- 
ing intensities up to 70%. The response curve of LAI versus grazing intensity was 
used in a mass balance model, called SNAP, of SOC dynamics based on previous 
data from the Serengeti. The model predicted SOC to increase at intermediate graz- 
ing intensity, but then to decline rapidly at the highest grazing intensities. The SNAP 
model predictions were compared with observed SOC stocks in the 24 grazed plots 
of a 10-year grazing exclosure experiment at eight sites across the park that varied 
in mean annual rainfall, soil texture, grazing intensity and plant lignin and cellulose. 
The model predicted current SOC stocks very well (-R 2 > 0.75), and suggests that 
compensatory plant responses to grazing are an important means of how herbivores 
might maintain or increase SOC in tropical grasslands. 
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INTRODUCTION 

Soil organic carbon (SOC) in grasslands and savannas represents one of the largest 
reservoirs of carbon on earth (Conant et ah, 2001; Lai et ah, 2007). Due to widespread 
unsustainable land use due to overgrazing, loss of native herbivores, excessive fires, and 
the potential to reverse the impacts of these uses, SOC is one of the more important 
potential sinks of greenhouse gases in the effort to mitigate climate change (Mannetje, 
2007; Ellis & Ramankutty, 2008; Smith et ah, 2008). A major question is what 
management practices in grasslands, in the form of grazing, fire, fertilization, 
re-vegetation and restoration, etc., can lead to net sequestration of carbon. 

Recent reviews demonstrate that herbivores can have dramatically different effects on 
soil organic carbon (SOC), both positive and negative, depending on soil type, 
precipitation, plant species composition and grazing intensity (Milchunas & Lauenroth, 
1993; Derner & Schuman, 2007; Pineiro et ah, 2010; McSherry & Ritchie, 2013). These 
results have been obtained largely in temperate grasslands grazed by livestock (McSherry 
& Ritchie, 2013), and relatively little is known about grazing impacts on SOC in tropical 
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grasslands. In the tropics and in certain temperate climates, plant species composition is 
dominated by warm-season (C4) grasses that, among other adaptations to grazing, invest 
heavily in rhizomes and other storage organs that allow them to respond quickly to both 
rainfall and defoliation (McNaughton, 1985; Milchunas & Lauenroth, 1993; Dubeux et ah, 
2007). Compensatory responses to grazing can involve sacrificing stems for leaves and 
thus maintain leaf area and potential carbon inputs despite carbon off-take by herbivores 
(Anderson, Dong & McNaughton, 2006; Zheng et ah, 2010; Ziter & MacDougall, 2013). 
Secondly, C4 grasses generally contain higher levels of lignin and cellulose (Barton et ah, 
1976), which are generally recalcitrant to decomposition, and produce extensive 
aboveground litter that frequently burns during extended dry seasons (Trollope, 1982; 
McNaughton, Stronach & Georgiadis, 1998; Sankaran, Ratnam & Hanan, 2008). Thirdly, 
benign temperatures allow for the prevalence of macro-decomposers, such as termites 
and dung beetles, to rapidly incorporate senesced plant material and herbivore dung into 
soil (Hanski & Cambefort, 1991; Freymann et ah, 2008; Dungait et ah, 2009; Risch, 
Anderson & Schutz, 2012). These features may dramatically alter herbivore impacts on 
SOC in tropical grasslands as compared to temperate grasslands with equivalent 
precipitation and soil type (McSherry & Ritchie, 2013). 

Of these three important features, the consequence of plant compensation for soil 
carbon is probably the least well understood. Leaf area index (leaf area/ area sampled), or 
LAI, has been shown to decrease in response to grazing in some contexts but not others 
(Anderson, Dong & McNaughton, 2006; Zhao et ah, 2009; Cruz et ah, 2010; Zheng et ah, 
2010; Ziter & MacDougall, 2013). Studies in Serengeti National Park (McNaughton, 1985; 
McNaughton, Milchunas & Frank, 1996) and west Africa (de Mazancourt, Loreau & 
Abbadie, 1999) have shown biomass compensation to grazing, but not compensation for 
lost leaf area specifically. Furthermore, the quantitative manner in which LAI changes as 
a response to varying grazing intensity has not been fully explored. A major question is 
whether LAI can be maintained or even increase with increasing grazing intensity and at 
what level of grazing intensity it declines. 

The response of LAI to herbivory may be a key driver of soil carbon response to 
grazing. Leaf area is highly correlated with gross photosynthesis and net carbon (C) 
assimilation (Craine et ah, 2001; Ellsworth et ah, 2004; Reich, Wright & Lusk, 2007). 
Maintenance of leaf area following defoliation potentially allows C inputs into important 
carbon pools, such as biomass and SOC, to be little affected by herbivory, while loss of 
leaf area can dramatically reduce C inputs. Because most current soil carbon dynamic 
models lack an explicit accounting of herbivore impacts on leaf area (Paustian, Parton & 
Persson, 1992; Schimel et ah, 1994; Parton et ah, 1995), the relative importance of leaf 
area-related C inputs relative to C losses from fire, decomposition, and herbivore 
consumption have not been explored. 

The goal of this paper is to show how leaf area index (LAI), or leaf area/area sampled, 
and plant allocation to stem and leaf changes with grazing intensity in a tropical savanna 
and then to explore how such possible compensatory responses might affect soil carbon 
dynamics. Hence, this paper (1) presents the results of a 2-year grazing exclosure 
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experiment in Serengeti National Park to measure effects of varying grazing intensity on 
plant LAI and proportion of biomass allocated to stems vs. leaves, (2) develops a model of 
soil carbon dynamics based on extensive past research in the Serengeti as a framework for 
evaluating the impact of LAI response to grazing, and (3) compares predictions of this 
model with observed soil carbon stocks across the eight sites of a 9-year grazing exclosure 
experiment in which key model input parameters were independently measured. 

The model of carbon dynamics, called SNAP (referring to its origin in Serengeti 
National Park) accounts explicitly for fixed carbon, as determined by LAI, soil texture 
and rainfall, and its allocation to roots versus shoots, and its fate: decomposition, 
consumption by herbivores, combustion in fire or residence in a long-term recalcitrant 
SOC pool. The model has five input variables: mean annual rainfall, grazing intensity, fire 
frequency, aboveground proportion of cellulose + lignin, and soil texture (percent sand), 
based on the suggested most important fates of stable carbon in the Serengeti ecosystem 
(McNaughton, 1985; McNaughton, Banyikwa & McNaughton, 1998; McNaughton, Stronach 
& Georgiadis, 1998; Holdo et al, 2007, Holdo et al, 2009). The model and its functional 
relationships are populated with parameter values available from a substantial literature 
on grazing, plants, soils, and microbes from Serengeti National Park plus some additional 
new data on grazing intensity, plant lignin and cellulose and soil moisture gathered from 
a long-term (9 year) grazing experiment across eight sites in Serengeti National Park that 
vary in grazing intensity, rainfall, soil texture, fire history, and plant species composition. 

To evaluate the role of plant compensation, a sensitivity analysis of model parameters, 
including grazing intensity, in the SNAP model is conducted. The model also is used to 
make deductive predictions of current soil carbon stocks under the assumption of a mass 
balance steady-state resulting from the persistence of mean rainfall, plant lignin and 
cellulose, soil texture, fire history, and grazing intensity conditions over an extended 
period of prior years. These predictions are then compared with observed soil carbon 
measured in grazed plots of the 9-year exclosure experiment using linear regression, with 
tests for overall model fit (R 2 ) and tests for differences in slope from a value of 1 and an 
intercept value of zero. This analysis provides an explicit and quantitative test of the 
hypothesis that plant compensatory responses are important in the influence of 
herbivores on soil carbon. 

This approach applies a deductive (Belovsky, 1984; Belovsky, 1994; Overmars, de Groot 
& Huigen, 2007) rather than inductive (Burnham & Anderson, 2004) evaluation of the 
impact of possible plant compensation. The use of a deductive model explicitly combines 
the logically-derived hypothesized mechanisms by which different factors might drive 
soil carbon dynamics based on prior knowledge of the Serengeti. A sensitivity analysis of 
the model assesses the relative importance of a given variable or parameter relative to 
others. The advantage of a deductive approach is that it avoids the uncertainty in 
inference derived from associating soil carbon dynamics with multiple, possibly 
auto-correlated, independent variables in an inductive approach (Burnham & Anderson, 
2004). Further, a deductive model builds in non-linearity of relationships, which is largely 
avoided in inductive approaches but may be likely to occur in the relationship between 
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LAI and grazing (McNaughton, 1985; de Mazancourt, Loreau & Abbadie, 1999; Anderson, 
Dong & McNaughton, 2006). 

MATERIALS AND METHODS 

Study area 

Most of the data presented in this study were collected in Serengeti National Park (SNP), 
Tanzania. SNP contains nearly 3 million migratory wildebeest (Connochaetes taurinus), 
zebra (Equus burchelli), and Thomson's gazelles (Gazella thomsoni), plus multiple species 
of resident herbivores > 10 kg in size, all distributed heterogeneously (Anderson et al, 
2010) over 25,000 km 2 (a mean density of 120/km 2 ) (McNaughton, 1985; Sinclair et al, 
2007). These herbivores consume a major fraction (mean 61.8% ± 12.2 SEM) of 
aboveground primary production (McNaughton, 1985). Vegetation is dominated by C4 
grasses, featuring Themeda triandra in more lightly grazed areas and Digitaria 
macroblephara and Pennisetum menzianum in more heavily grazed areas. The Serengeti 
region features a pronounced rainfall gradient driven mixed somewhat independently 
with a soil gradient of silty, organic matter-rich soils in the southeast to clay soils in the 
west and sandy, low nutrient soils in the north (Anderson et al, 2007; Sinclair et al, 2007; 
Anderson et al, 2010). The clustered distribution of resident herbivores (Anderson et al, 
2010) and variable movements of migratory species results in highly variable historical 
grazing intensity across the landscape. Strong preferences by herbivores for plants found 
at the top rather than bottom of slopes plus a tendency to avoid the edges of woodlands 
results in variable grazing intensity over distances of <50 m (Gwynne & Bell, 1968; 
Anderson, Dong & McNaughton, 2006; Anderson et al, 2007). 

Design 

Plant compensation to grazing was assessed near the Serengeti Wildlife Research Centre 
near Seronera (34° 50' E and 2° 25' S) inside the Park. The area is open Acacia tortilis 
savanna woodland, with dominant grasses Themeda triandra and Digitaria 
macroblephara, and hosts resident herds of impala (Aepyceros melampus), buffalo 
(Syncerus coffer), and Thomson's gazelles (Gazella thomsoni), among others, and is 
routinely visited by the annual wildebeest and zebra migration. Plant responses to grazing 
were measured inside and outside of grazing exclosures established 2 years previously in a 
10 ha open grassland with less than 5% tree cover near the Serengeti Wildlife Research 
Centre (2-year exclosure experiment). The proportion of biomass as stem versus leaf was 
measured for all plants clipped from six randomly selected 15x15 cm quadrats, three 
inside and three outside each of twelve exclosures of various size (10-200 m 2 ) erected 
between 2009 and 2010 to protect experimental gardens or nutrient addition 
experiments. LAI, proportion of stem, and aboveground biomass were measured in the 
unused "buffer zone" just inside each fence and then in quadrats 4 m perpendicular to 
each fence paired with each quadrat inside the fence. The objective was to encounter 
different localized grazing intensities that resulted from sampling beneath thorny shrubs 
or other plant protections versus patches with highly preferred forage species. 
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To generate additional input data for a model of soil carbon dynamics, a longer-term (9 
year) grazing exclosure experiment established by Sam McNaughton in 1999 (Anderson, 
Ritchie & McNaughton, 2007) was used to provide data on soil moisture as a function of 
rainfall, and independent measures of grazing intensity, plant lignin and cellulose, and 
SOC to test the SOC dynamic model. A line of six plots (4x4 m), spaced 10 m apart in a 
line transect were established in October 1999 at each of 8 sites located 10 or more km 
apart, in grassland areas visited almost entirely by grazing, rather than browsing, ungulate 
species (Anderson, Ritchie & McNaughton, 2007). Sites were chosen to be within 1 km of 
large, permanent concentrations of grazing herbivores (Anderson et al, 2010) and varied 
considerably in annual rainfall, soil type, and fire frequency Three randomly selected 
plots within the transect at each site were fenced with 2 m high, 8 cm mesh wire. These 
fences effectively excluded all grazing mammals > 10 kg, since animals preferred to go 
around rather than jump fences. The remaining three plots at each site were unfenced 
controls. Woody vegetation at our experimental sites was typically sparse (less than 10 
woody plants/ha) but herbaceous legumes, such as Indigofera volkensii, were ubiquitous 
and reasonably abundant ( > 2% cover) at all our study sites. Our eight sites included two 
on the treeless shortgrass plains of southeastern SNP, and the other six sites were 
representative of woody savannas broadly distributed across SNP (Ruess & Seagle, 1994). 

Measurements 

The effect of grazing on plant (grass) allocation to leaf versus stem and on leaf area index, 
LAI (cm 2 leaf/ cm 2 ground area) was assessed. Clipped plants for all plant species in a 
quadrat were taken back to the lab and, while fresh, laid out and traced on graph paper, 
and then leaves were clipped from stems and stored, dried at 45°C, and weighed 
separately (as with field biomass). The proportion of leaf for each quadrat was calculated, 
along with LAI (total leaf area (cm 2 ) divided by quadrat area (cm 2 )). 

Another key variable, GI (grazing intensity) varies locally and was estimated in both 
the 2-year and 9-year experiments as the fractional difference in standing aboveground 
biomass (AGB) : GI = 1 - AGB^/ABG^ (Anderson, Ritchie & McNaughton, 2007). GI is 
an indirect measure of the fraction of annual production consumed by grazers 
(McNaughton, 1985; McNaughton, Milchunas & Frank, 1996). AGB was measured by 
clipping and weighing all aboveground plant material, excluding litter, after drying at 
45°C for three days to constant mass. AGB was measured in the 2-year exclosure 
experiment in the paired quadrats inside and outside exclosures used to measure LAI. 
The quadrat inside the exclosure was assigned GI = 0, and the paired quadrat outside the 
exclosure was assigned GI based on the difference in AGB between the paired quadrats. 
In the 9-year exclosure experiment, biomass was clipped from four 15 x 15 cm quadrats 
within each grazed and exclosure plot in late May or June at the end of the wet season in 
2006, 2009, and 2010. Biomass measured in each of the three grazed plots at each site was 
compared to biomass in each of their nearest neighbor exclosure plots (Anderson, Ritchie 
& McNaughton, 2007), yielding three estimates of grazing intensity at each site. 

To measure lignin pluscellulose fraction (LIGCELL), dried plant material from the 
entire clipped quadrat was ground through a 0.9 mm (40 mesh) screen in a Wiley mill 
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and then subjected to a sulfuric acid hydrolysis method (Sluiter et al, 2010), i.e., 
sequential digestion in first neutral detergent manufactured by Ankom Technology Corp., 
Macedon, New York, USA (sodium lauryl sulfate, 50.02%, EDTA disodium, dehydrate, 
31.03%., sodium borate, decahydrate 11.35%, and sodium phosphate, dibasic, anhydrous, 
7.6%), followed by Ankom manufactured acid detergent (20 g cetyl trimethylammonium 
bromide (CTAB) diluted in 1 L 1.00 NH 2 S0 4 ), and then in 72% sulfuric acid, followed 
finally by ashing in a muffle furnace at 500° C for 24 h. Lignin and cellulose fractions were 
inferred by the mass lost between the remainder following neutral detergent digestion 
and the remainder after ashing, expressed as a proportion of the original sample mass. 

Annual rainfall over the 9 years of the study was obtained from rain gauge data at 23 
sites collected by the Serengeti Ecological Monitoring Program (Anderson et al, 2007). 
Fire frequency was determined by interpreting MODIS satellite visible and near infra-red 
imagery (Dempewolf et al, 2007) for July-October for 2000-2008 to generate maps of 
fire occurrence for each year with a resolution of 250 m. Grazing experiment site GPS 
locations were overlaid with fire maps for each year in a GIS (ArcGIS 9.1 byESRI®) to 
determine the number of times a site burned over the period 2000-2008. 

Effects of precipitation on soil moisture and number of days in which microbes were 
likely to be active (WETDAYS) were determined from field measurements of soil 
moisture in the long-term experiment (unpublished data generously provided by Sam 
McNaughton). Five cm diameter soil cores were taken to a depth of 10 cm each month 
from December 1999 to June 2002 (a total of 31 days) at each of three grazed plots at the 
eight sites of the long-term Serengeti grazing experiment. Moisture content was 
determined gravimetrically and monthly rainfall was measured from rain gauges 
mounted on 1 m tall fence posts at each site. A recent review (Manzoni, Schimel & 
Porporato, 2012) suggests that, in semi-arid soils, microbial activity ceases when water 
potential drops below — 14 cm H2O. Water potential was converted into soil moisture, 
calibrated for the texture of Serengeti soils (Table 1) (Cosby et al, 1984) and suggested 
that Manzoni, Schimel & Porporato (2012) threshold corresponded on average to 10% 
gravimetric water content. The number of days at each of the eight sites that soil moisture 
was 10% or higher was determined and then used to build a regression to estimate the 
number of days per year with active soil microbes as a function of annual rainfall. 

As a test of the new soil carbon model, SOC was measured at the same time as biomass 
in each plot from 8.3 cm diameter x 40 cm deep cores from each of the three biomass 
quadrats per plot. Each core was spread on a portable table top and crushed with a rolling 
pin to break up soil aggregates. All pebbles, rocks, and visible roots were removed, and 
then soil was sieved through 0.7 mm mesh. Cores of known volume were pooled, stored 
in sealed plastic bags, and then dried for 6 days at 45°C and then weighed to determine 
bulk density (g soil/ cm 3 ). Both SOC and nitrogen contents were determined with the 
Walkley- Black method at the Soils Analysis Laboratory at Sokoine University of 
Agriculture in Morogoro, Tanzania. Reported error was ±0.004 g/gC. Soil composition 
(sand, silt and clay fractions) was determined by a micropipette method. Soil carbon 
density was measured by multiplying SOC content by bulk density (g/m 2 ) to 40 cm depth 
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Table 1 Variables and key parameters in the SNAP model, their units, and sources of data. 



Parameter 


Units 


Type 


Source 


ANPP max 


— =3 — =n — 

CI TT"I ~\TY 

g m yr 


T 1 t~£*V3 t~l 1 f£» 

-LillCI dLUIC 


IVlctAllllUlll dUUVCl'lUUllU UIUUUCLUJll ylvlLl\lilAvillUfl } lyOD) 


ANPP est 


—2 -1 
g m yr 


IVHJU.C1 


(4 to 7pr_ tyi /"\ci i i i A \TPP 
Ul aLCL illULLlLiCLl ril\rr 




—2 —1 
gm yr 


Literature 


Belowground production {McNaughton, Butiyikwa & NLc- 

IvtluVrllurl, lyyo) \r\Q. jd) 




/- —2 -1 
g(-m F 


iviooei 


Dung-derived SOC 


FIRE 


#/yr 






GI 


proportion 


Measured 


1- (grazed biomass/ungrazed biomass) (Table 4) 


LAI 




Measured 


Leaf Area Index (dimensionless); This study (Fig. 2B) 


LIGCELL 


proportion 


Measured 


Cellulose + lignin (%) aboveground biomass, (Table 4) 


MRESP 


gCm" 2 d" 1 


Literature 


(Ruess & Seagle, 1994) (Fig. 3D) 


P L 


proportion 


Measured 


Proportion leaf; This study (Fig. 2A) 


PDSOC 


-2 -1 
gC m yr 


Model 


Plant-derived SOC 


RAIN 


mm/yr 


Interpolated 


Mean annual rainfall (Anderson et al„ 2010) (Table 4) 


SAND% 


percent 


Measured 


From a standard soil analysis, measured (Table 4) 


WETDAYS 


dyr" 1 


Measured 


#days with soil > 10% water; This study (Fig. 3C) 


WHC 


proportion 


Literature 


Water holding capacity (Ruess & Seagle, 1994) (Fig. 3A) 



Table 2 Parameters for Monte Carlo simulations. Parameters and their errors used in predicting soil carbon 
stocks in Serengeti grasslands and in the Monte Carlo simulations. 

Parameter Input Coefficient 







Bo 


SE 


Bi 


SE 


B 2 


SE 




GI 


0.597 


0.061 


0.24 


0.023 






LAI 


GI 


1.15 


0.027 


-0.015 


0.0011 


4.6 


0.2 


ANPP max 


RAIN 


-27.5 


17.6 


0.84 


0.07 






WHC 


SAND% 


1.04 


0.05 


-0.0070 


0.0008 






BNPP est 


RAIN 


958.8 


165.0 


-0.82 


0.27 






WETDAYS 


RAIN 


-0.025 


0.033 


0.00043 


0.00006 






MRESP 


SOC 


-0.58 


0.45 


0.00044 


0.00007 






CARBON 


N/A 


0.45 


0.02 











of the pooled cores in each plot. A 40 cm depth was used, as this typically corresponds to 
the depth of >90% of roots (McNaughton, Stronach & Georgiadis, 1998), and multiple 
sites featured a hardpan or dense clay layer at 40-50 cm that prevented reasonable 
sampling to greater depths. 

Literature data and re-analyses 

Literature data, mostly from prior studies in Serengeti, were re-analyzed to specifically 
estimate soil water holding capacity (WHC), aboveground production, belowground 
production, and soil microbial respiration of carbon. For sources and units of literature 
data, see Table 1. The relationships used are reported in the Results and Discussion 
(Table 2). 
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Figure 1 Hypothetical major fates of carbon in tropical grassland, as the basis for a practical soil 
carbon dynamic model. Net fixed carbon becomes resident soil organic carbon (SOC) through two ma- 
jor pathways: plant-derived SOC not consumed by microbes (in soil or the guts of macro-decomposers 
like termites) (heavy solid arrows), and dung-derived SOC (heavy short-dashed arrows) not assimilated 
by grazers and gut or soil microbes. All other carbon is combusted through fire (black arrow) , or respired 
by grazers and microbes in soil or the guts of macro-decomposers (long dashed arrows). 



Model development 

The flow of fixed carbon through hypothetical tropical grassland was described according 
to four criteria. First, mass balance was assumed, so that all fixed carbon resided in pools, 
such as biomass or SOC, or suffered the fate of being combusted, consumed, or respired. 
Second, this flow was described in the minimum number of fluxes and pools necessary to 
account for the fate of carbon (Fig. 1). Third, the most parsimonious set of input variables 
that could describe the changes in these pools was used. Finally, a time step of one year 
with a distinct wet and dry season each year was assumed, with production and grazing 
resulting in a standing biomass at the beginning of the dry season each year that 
subsequently either combusted if there was fire or was otherwise decomposed. More 
complicated time dynamics of daily or seasonal changes have been used recently (Holdo 
et at, 2007, 2009), but this introduces possibly unnecessary complexity. Above- and 
belowground production are therefore seasonal accumulations of biomass as modified by 
average grazing intensity over the course of the wet season. Biomass not consumed by 
herbivores may then burn at some frequency during the dry season. These assumptions 
led to a relatively simple set of pathways for stable carbon (Fig. 1): (1) input into soil 
organic matter from decomposing shoots or roots, (2) combustion in fire, or (3) input 
into soil organic matter through dung. All non-lignin or cellulose compounds were 
assumed to turnover within 1-2 years or less (Ruess & Seagle, 1994) and were 
consequently ignored. 

The model was based on four major assumptions. First, tropical grasses potentially 
exhibit compensatory responses to defoliation that maintain similar leaf area, and thus 
photosynthetic capacity, across a broad range of grazing intensities (McNaughton, 1985). 
Second, the largest carbon inputs to the soil organic matter pool occur through 
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decomposition of above- and belowground biomass and through incorporation of 
herbivore dung into soil. Third, the major losses of SOC derive from combustion (fire), 
herbivore respiration and soil microbial respiration. Finally, all plant C in plant tissue 
other than lignin and cellulose was assumed to be assimilated and respired by herbivores, 
microbes in soil or the guts of macro-decomposers. Lignin and cellulose are the most 
recalcitrant forms of carbon to decomposition, and likely account for the majority 
(>90%) of long-lived SOC (Ganjegunte et ah, 2005; Frank, Pontes & McFarlane, 2012). 

The SOC (g/m 2 to 40 cm depth) measured in 2009 in the grazed plots in the long-term 
grazing experiment represented SOC e q under the reasonable assumption that conditions 
in these plots had been approximately similar in terms of wildlife use, plant species 
composition, and climate since 1979 {Sinclair et ah, 2007; Holdo et ah, 2009). The average 
of grazing intensities measured in 2006, 2009, and 2010 was assumed to reflect historical 
long-term grazing. The other input variables either were already a long-term average 
(RAIN, FIRE) or are unchanged by land use (SAND%). 

Statistics 

The association between LAI and grazing intensity (GI) was tested with linear regression 
in SPSS version 19 (IBM, Armonk, New York, USA). Non-linear relationships were also 
explored, and based on the appearance of the data, the R script "nls" was used with a 
function Bo — Bi exp(— B2 GI), where Bj are fitted coefficients. Linear regression in SPSS 
was used to build regression relationships from literature data (Table 2) and to compare 
observed SOC to predicted SOC eq from the SNAP model. In the comparison of the 
observed and predicted SOC, the hypotheses that the intercept of the regression line was 
equal to zero and that the slope was equal to one were evaluated with z-tests. 

The SNAP soil carbon dynamic model was constructed from a number of empirical 
relationships that contain parameters (intercepts, slopes, etc.) estimated with error. A 
Monte Carlo analysis {Ogle et ah, 2010) was conducted by sampling the normal 
distributions inferred from measured model input parameters and parameters from 
relationships used to construct the model. This process was repeated 100 times to 
generate a mean and s.e.m. for the estimated SOC eq for each site and each paired grazed 
and ungrazed plot in the 9-year exclosure experiment. 

RESULTS AND DISCUSSION 
New and prior data analyses 

The proportion of biomass in leaves increased significantly (P < 0.001) with increasing 
grazing intensity (Fig. 2A). In concordance with this pattern, Leaf Area Index (LAI) 
declined significantly with increasing grazing intensity (GI) according to linear regression 
(R 2 = 0.36, df = 1,34, P = 0.01). However a non-linear function, LAI = B 0 - B x exp 
(B 2 GI) fit the data much better (R 2 = 0.58, df = 1,33, P < 0.001) and visibly matched 
better the consistently high LAI found at low to intermediate grazing intensity and the 
dramatic decline in LAI above a grazing intensity of 70% (Fig. 2B). These two patterns 
support the hypothesis that plants maintain C inputs to the Serengeti ecosystem despite 
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Figure 2 Plant compensatory responses to grazing in Serengeti. Key relationships for assessing plant 
compensation to grazing: (A) proportion of leaves (Pi) and (B) Leaf Area Index (LAI, cm 2 /cm 2 ) of 
grazed and ungrazed plants, each as a function of grazing intensity (GI). 



heavy removals of biomass by grazers. The maintenance of LAI across such a broad range 
of grazing intensities contrasts with the typical pattern of dramatic reductions in LAI seen 
mostly in temperate systems under relatively moderate grazing intensities (30%-60%) 
(Zhao et ah, 2009; Cruz et al, 2010; Zheng et al, 2010). The pattern corroborates results 
of enhanced productivity under grazing (McNaughton, 1985; Milchunas & Lauenroth, 
1993; Frank, Kuns & Guido, 2002) that are at least partly derived from compensatory 
morphological changes by plants to grazing (Anderson, Dong & McNaughton, 2006; Ziter 
& MacDougall, 2013). Despite these previous results, it is still not clear why the C4 grasses 
typically dominant in the Serengeti and elsewhere respond differently to moderately 
intense grazing than C3 grasses (Derner & Schuman, 2007; McSherry & Ritchie, 2013). 

As expected, water holding capacity exhibited a strong, significant (P < 0.001) 
negative association with increasing sand content (Ruess & Seagle, 1994) (Fig. 3A), and 
root production showed a significant (P = 0.0015) decline in association with higher 
rainfall (McNaughton, Banyikwa & McNaughton, 1998) (Fig. 3B). The proportion of days 
soil held more than 10% moisture increased significantly (P = 0.006) with mean annual 
rainfall across the eight long-term grazing experimental sites (Fig. 3C). Finally, maximum 
microbial respiration rates were strongly and significantly (P < 0.001) associated with 
existing SOC stocks across the Serengeti landscape (Ruess & Seagle, 1994) (Fig. 3D). 
Regression equations are reported in the figures and were used in the development of the 
SOC dynamic model. 

SOC dynamic model 

Based on the results for LAI response to grazing (Fig. 2B), the SOC model for Serengeti 
National Park (hereafter referred to as the SNAP model) was constructed with five input 
variables: RAIN, mean annual rainfall (mm/yr); GI, grazing intensity (l-(grazed 
biomass/ungrazed biomass)) (McNaughton, 1985); FIRE, fire frequency or the number of 
fires recorded, divided by the number of years over which fires were monitored; 
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Figure 3 Functional relationships used to build the SNAP soil carbon model. Key analyses of liter- 
ature data to build empirical relationships for key model parameters. (A) the decline in water-holding 
capacity of soils with increasing sand content (Ruess & Seagle, 1994), (B) annual belowground produc- 
tion as a function of rainfall (McNaughton, Stronach & Georgiadis, 1998), (C) number of moist soil days 
(WETDAYS) as a function of annual rainfall (RAIN), and maximum microbial maximum respiration 
rates as a function soil carbon stocks (SOC) (from Ruess & Seagle (1994)). 



LIGCELL, lignin and cellulose content (proportion) in aboveground herbaceous 
vegetation, and SAND%, (McNaughton, 1985). The sources of these data are summarized 
in Table 1 and specific regression equations with s.e.m. error terms are provided 
in Table 2. 

Two additional factors that can be important influences on carbon dynamics, nitrogen 
(Pineiro et ah, 2010) and soil temperature (Miller et ah, 2004), were also considered but 
not included as inputs into the SNAP model. Production in the Serengeti is assumed to be 
co-limited by nitrogen and water, and total soil N was included in a validated model of 
plant productivity for Serengeti (Holdo et ah, 2007). Simplifying the functional forms 
used in Holdo et al. s model and solving at steady state revealed that ANPP max is likely to 
be proportional to the product of soil N and rainfall. However a multiple regression of the 
logarithms of ungrazed biomass (production) in the long-term grazing exclosures against 
the logarithms of soil N (ln(SoilN) and mean annual rainfall (ln(RAIN)) revealed that 
ln(RAIN) was a highly significant factor, but ln(Soil N) was not (P = 0.45). This outcome 
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suggests that N is a far less important factor influencing productivity in the Serengeti 
than is rainfall. In addition, soil N or soil C:N was not a significant variable explaining 
residual variation in microbial respiration in Serengeti soils following regression if SOC 
was included in a regression model (partial R 2 = 0.05, N = 16, P < 0.41). Finally, the 
influence of other variables that are often correlated with higher N, such as finer-textured 
soils (Ruess & Seagle, 1994) and more intense grazing (Holdo et at, 2007) will already be 
included in the model in other relationships. 

The influence of soil temperature on microbial respiration was also not incorporated 
into the model. Temperature in most tropical soils is not as variable as in temperate 
systems, typically ranging over only 3°C over an entire year (Caquet et ah, 2012), and 
increased soil temperature is typically accompanied by drier conditions, thus canceling 
any temperature-related increases. Consequently, temperature was viewed as a 
difficult-to-measure parameter that would likely add relatively little explanatory power to 
the model. 

Fixed carbon was accounted for as biomass produced during the wet season, and was 
assumed to be allocated either aboveground or belowground. Estimated aboveground 
production in year r, ANPPf \ was determined as potential production, modified by 
grazing impacts on proportion biomass as leaf and specific leaf area. Potential production 
as a function of rainfall (Table 2) was derived from McNaughton (1985) empirical 
relationship and weighted by the adjusted water holding capacity (AWHC), as 
determined by soil sand content AWHC = 1.33 - 0.0075*SAND%. This is the 
relationship between water holding capacity (WHC) and SAND% derived from Ruess & 
Seagle (1994) (1.08 - 0.0075*SAND%) but with 2 standard deviations in their measured 
WHC added to the intercept. This converted WHC into a coefficient that adjusted 
McNaughton (1985) productivity relationships, which reflected average soil properties, for 
whether soil is higher or lower than average in sand content. Adjusting potential 
productivity for water holding capacity incorporates how soil influences soil moisture for 
plant growth at a given rainfall amount. 

Thus: 

ANpp max _ (Q 84 # RAm _ 27 5) # (1 3 3 _ 0i0 075 * SAND%). (1) 

Potential aboveground production is then adjusted by leaf area index, LAI, to calculate 
estimated productivity ANPPf \ which accounts for the stimulation of leaf area at the 
expense of loss of stem due to grazing. 

ANPPf 1 = LAI * ANPP™ ax . (2) 

As measured by McNaughton (1985), ANPPf 1 can be higher than ANPPf at low to 
intermediate grazing intensities because grazing keeps plants at a faster-growing stage 
during the wet season. At high grazing intensities, leaf area declines and ANPPf 1 
approaches zero. LAI was determined from the non-linear empirical relationship for LAI 
measured in this study (Fig. 2B): 

LAI = P L /0.6 - 0.015 * EXP(4.6 * GI), (3) 
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where Pl is the proportion of biomass as leaves, GI is grazing intensity, and EXP is the 
exponential function. Pl was also estimated from plant compensation to grazing 
measurements (Fig. 2A) 

P L = 0.6 + 0.24 * GI. (4) 

The ratio Pl/0.6 in Eq. (3) incorporates the grazing-induced adjustment in the 
proportion of leaves into the calculation of LAI, such that LAI = 1 at GI = 0 and 
Pl = 0.6, the y- intercept value in Fig. 2A. Combining Eq. (2)-(4) yields estimated 
aboveground productivity as a function of rainfall, soil texture, and grazing intensity: 

ANPP t est - {[(0.84* RAIN -27.5)* (1.33 - 0.0075 *SAND%)] 

* [(0.6 + 0.24 * GI)/0.6 - 0.015 * EXP(4.6 * GI)]. (5) 

Belowground production in year t generally declines with increasing rainfall 

(McNaughton, Banyikwa & McNaughton, 1998) (Fig. 1C) 

BNPP t est = 917.4-0.763* RAIN. (6) 

With these C input functions, inputs to SOC were calculated from the two major 
sources, plant-derived carbon PDSOC t (solid arrows in Fig. 1) and dung-derived carbon, 
DDSOC t (thick short-dashed arrows in Fig. 1). Plant-derived carbon consists of the 
carbon fraction of lignin and cellulose in roots plus the carbon fraction of lignin and 
cellulose in aboveground biomass that is not consumed by grazers or burned. Total 
non-ash-free carbon content of aboveground plant tissue (Tao et ah, 2012) averages 0.45. 
Lignin and cellulose in roots was not measured but was estimated to be 5% higher than 
that in aboveground tissues (Andrioli & Distel, 2008; Semmartin, Garibaldi & Chaneton, 
2008). Grazing and fire interact in that fire can only burn biomass not consumed by 
grazers, so that aboveground decomposition is only of the fraction of biomass not grazed 
(1-GI) and not burned (1-FIRE). Fire and grazing were assumed to have no important 
effects on belowground production (McNaughton, Banyikwa & McNaughton, 1998). 

PDSOCt = 0.45 * [LIGCELL * ANPP t est * (1-GI)* (1-FIRE) 

+ (LIGCELL + 0.05) * BNPP t est ] . (7) 

Dung-derived carbon was assumed to be the lignin and cellulose fraction of biomass 
consumed by grazers. All stable C in dung was assumed to be incorporated in soil because 
of the activities of termites and dung beetles, which often incorporate dung into soil 
within 2-5 days (Freymann et ah, 2008; Risch, Anderson & Schutz, 2012), allowing little C 
to volatilize. While a significant fraction of the dung is consumed and respired by termites 
and dung beetles, the stable C fraction likely will remain intact as soil organic matter and 
thus contribute to SOC. As for aboveground plant tissue, C content of dung was 0.45. 

DDSOCt - LIGCELL * 0.45 * GI * ANPP t est . (8) 
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Carbon losses associated with microbial respiration (MRESP t ) were accounted for as a 
function of the number of days in which soil moisture exceeds a threshold for microbial 
activity, which was determined to be 10% (see Methods). The proportion of days during 
the typical 240 day wet season in which microbes would be active (WETDAYS) was 
estimated from the average annual rain fall from the regression in Fig. 3C: 

WETDAYS = (0.00044 * RAIN - 0.025) * 240 (9) 

WETDAYS can be multiplied by the maximum rate of microbial respiration (Fig. 3D), as 
determined by lab incubations of Serengeti soils: 

MRESPt = WETDAYS * (0.7 + 0.3 * SAND%/100) * (0.00044 * SOQ - 0.579), (10) 

where the expression (0.7 + 0.3 * SAND%/100) corrects the respiration rate for greater 
accessibility to SOC in sandy soils, as suggested by data in Ruess & Seagle (1994). 

These equations account for all the major inputs and losses of stable C into soil. The 
change in sequestered carbon ( ASOCt) can be estimated by adding the plant derived SOC 
(PDSOQ) to the dung derived SOC (DDSOQ) and then subtracting the carbon lost 
through microbial maintenance respiration (MRESPt) as follows: 

ASOCt = PDSOCt + DDSOCt - MRESPt. (11) 

By setting ASOCt = 0, the above equation can be solved for the SOCt term in MRESPt to 
yield an equilibrium SOC eq . 

SOCe q = [PDSOCt + DDSOCt + WETDAYS * (0.579) 

* (0.7 + 0.3 * SAND%/100)]/[(0.00044 * WETDAYS 

* (0.7 + 0.3 * SAND%/100)]. (12) 

Model predictions 

The equation for SOC eq allows for a relatively straightforward analysis of how different 
factors influence soil carbon dynamics. It can predict ASOCt from knowledge of some 
initial soil carbon stock Eq. (11), or it can predict the equilibrium or steady-state carbon 
stocks following an extended period of certain conditions (grazing, fire, plant species 
composition, etc.) Eq. (12). My focus here is to explore the influence of different factors 
on equilibrium carbon stocks, as these represent different states that grassland systems 
might be approaching from current or past states. This analysis will also reveal the 
qualitative impacts of different factors on SOC stocks or fluxes. 

In the SNAP model, SOC eq increased with higher plant lignin and cellulose content, 
but decreased with increasing rainfall and soil sand content (Fig. 4). Increasing rainfall 
decreased C inputs through roots Eq. (6)) and increasing sand reduced water holding 
capacity and aboveground production as well as increased soil microbial respiration 
through increased access to organic carbon in coarser soils. In contrast to these 
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Figure 4 Predictions of soil carbon stocks in the Serengeti. Model predictions of SOCeq stocks to 40 
cm depth as a function of the key input variables. These include grazing intensity and fire frequency for 
different combinations (A-D) of low (RAIN = 450) and high (RAIN = 800) rainfall and low (SAND% 
= 25, fine-textured soils) and high (SAND% = 65, coarse-textured) soil sand content. Lignin and cellu- 
lose (LIGCELL) was set at 0.3 following Barton et al. (1976). 



straightforward effects, grazing impacts on SOC eq depended on fire frequency (Fig. 4). At 
low fire frequency, grazing reduced standing biomass and C inputs through plant 
decomposition, so SOC e q declined very little at lower grazing intensities but exhibited 
large decreases in SOC e q at grazing intensities >70%. At high fire frequency, high grazing 
intensities produced the same dramatic drop in SOC, but low grazing intensity led to 
lower SOCeq as compared to intermediate grazing intensities. 

Compensatory maintenance of LAI across a wide range of grazing intensities (Fig. 2B) 
featured prominently in these predictions because compensatory growth was predicted to 
maintain C inputs and aboveground biomass fuel for fires even when grazers consume 
the majority of available biomass. Other work in Serengeti shows that the shift from 
stems to leaves is chronic as the result from long-term replacements of taller more erect 
species with shorter, more leafy species (McNaughton, 1985; Anderson et al, 2007). As 
several papers outside Serengeti have shown, C assimilation, which is tightly coupled to 
leaf area, appears to be the major driver of carbon inputs into ecosystems (Powell et al, 
2006; Jarlan et al, 2008; Cleverly et al, 2013). Biomass reflects the ratio of C assimilation 
to turnover, but turnover can reflect allocation from shoots to roots, exudation of C to 
soil, consumption by herbivores, etc., and biomass is thus not directly proportional to C 
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Table 3 Parameter values used in the sensitivity analysis. 



Level 


Rainfall 


Fire frequency 


Grazing intensity 


Lignin and cellulose 


% Sand 


Low 


450 


0.25 


0.25 


0.15 


25 


Medium 


650 


0.50 


0.60 


0.25 


45 


High 


900 


0.75 


0.90 


0.35 


70 




LIGCELL SAND% 



Figure 5 Sensitivity analysis for the SNAP soil carbon model. Model sensitivity analysis of the mean 
(+SE) absolute value of percent change in model-predicted equilibrium SOC (SOCeq). Change results 
from variation in each input parameter by +10% while holding other parameters constant. Sensitivities 
(% change in equilibrium SOC ± SE) were calculated for low, medium and high initial values of each 
parameter (Table 3). 



inputs into the ecosystem and especially not the soil. Consequently LAI compensation 
following herbivory can be reasonably associated with sustained C inputs. 

Sensitivity and uncertainty analysis 

The sensitivity of calculated carbon stocks (% change) to 10% changes in parameter 
values (Table 3) shows that the model predictions are most sensitive to changes in grazing 
intensity especially at the upper range of grazing intensities (Fig. 5). ASOQ was about 
equally sensitive to changes in rainfall, fire frequency and lignin and cellulose content, 
and was least sensitive to changes in soil sand content. 

The Monte Carlo analysis yielded consistently small standard errors in estimated 
SOCeq of around 250-400 g/m 2 (around 4%) and an overall uncertainty of around 
14% (Fig. 6). These figures are consistent with the sensitivity analysis (Fig. 5) and 
suggest that errors in SOC predicted by the model are much less than the range of 
SOC (3000-9000 g/m 2 ) observed across the range of conditions in the Serengeti. 

Model validation 

Predicted SOC eq was compared to observed SOC (SOC 0 b s ) in two ways. First, the average 
predicted SOC eq and SOC 0 b s were compared at the site level (N = 8), i.e., averaged over 
all three plots at a site, which reflects that RAIN, FIRE, and SAND% were the same for all 
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three plots. Secondly, to test whether accounting for local variation in grazing intensity 
and lignin and cellulose (Table 4) affected model success, SOC e q and SOC G bs were 
compared at the plot level (N = 24). At the site level, average SOC eq at a site strongly 
predicted the mean SOC among grazed plots at that site (Fig. 6A): SOC 0 bs = 0.817 
(±0.12 s.e.m.) SOC eq + 858.1 (±1178.3 s.e.m.). The slope was not significantly different 
from 1 and the intercept was not different from zero (f -tests, P > 0.5). Likewise, at the 
plot level, SOCeq for a plot strongly predicted the mean SOC for that plot (Fig. 6B): 
SOCobs = 0.77 (±0.11 s.e.m.) SOC eq + 1122.1 (±1101.3 s.e.m.). However, even though 
the intercept was not different from zero the slope was significantly different from 1 and 
(z-tests, P > 0.5). Both slopes were less than 1, which suggests that the model may 
over-predict SOC under conditions that sustain the highest SOC (drier conditions with 
finer-textured soils, few fires, and moderate grazing intensities). Very few of the microbial 
respiration rates measured by Ruess & Seagle (1994) occurred under these conditions, and 
it is possible that respiration rates were higher than those incorporated in the SNAP 
model, which would result in lower than predicted SOC at equilibrium. Further 
exploration of plant compensatory responses, soil microbial respiration and other key 
processes under drier conditions is likely warranted. 

The strong correlation with observed carbon stocks (Fig. 6) suggests that the model 
captures the most important factors influencing soil carbon for the Serengeti system 
(Fig. 2) and that plant compensatory responses are a key component of soil carbon 
dynamics in the Serengeti ecosystem. The model is balanced in its sensitivity to different 
input parameters (Fig. 5), with the exception of its sensitivity to grazing intensity above 
70% (Table 3, Figs. 4 and 5). The model has enhanced sensitivity to grazing above this 
level because LAI declines rapidly under increasing grazing intensity once grazing 
off-take exceeds plants' ability to trade stem (Fig. 2B) or rhizome (not measured here) to 
maintain leaf area and hence carbon inputs to the ecosystem. 

As is implied by the solution for SOC e q (12) soil carbon stocks are predicted to increase 
linearly with lower soil sand content, lower rainfall, and greater plant lignin and cellulose. 
The predicted response of SOC e q to decrease with increasing rainfall appears 
counterintuitive but results from the decline in root biomass and productivity with 
increasing rainfall that was observed in the one available tropical grassland study 
(McNaughton, Banyikwa & McNaughton, 1998). This may occur because, within tropical 
grasslands, at even the lowest rainfall amounts (400-500 mm), aboveground productivity 
and biomass in the absence of grazing can exceed 600 g/m 2 (Table 4) enough to induce 
light limitation at the soil surface (Anderson et al, 2007). Consequently, plants may 
potentially trade off allocation to root production to compete for light aboveground. This 
may lead to declining root-derived carbon inputs at higher rainfall (Fig. 3B), and, since 
roots of herbaceous plants generally contain higher lignin contents than their 
aboveground tissues, reduced overall stable carbon inputs to soil. 

A key and unique aspect of the SNAP model is the interaction between fire and grazing 
(Fig. 4), which leads to a hypothesized peak in possible carbon stocks at intermediate 
grazing intensities (Fig. 4) that is amplified by more frequent fires. This is not a dynamic 
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Figure 6 Validation of the SNAP soil carbon model Model-predicted equilibrium soil organic carbon 
stocks (g/m 2 to 40 cm depth) compared to (A) mean (is.e.m.) observed soil carbon stocks (N = 3) 
measured at the eight 9-year grazing experiment sites and (B) the 24 grazed plots in the 9-year graz- 
ing experiment. Horizontal error bars were derived from Monte Carlo simulations based on parameter 
errors presented in Table 3. 



that is captured explicitly by previous models (Paustian, Parton & Persson, 1992; Parton 
et al, 1995; Pineiro et al, 2010), and is captured here specifically because of the influence 
of plant compensation. As expected, at high grazing intensities, the loss of potential fuel 
diminishes the impact of fire on SOC in the model. 

Literature reviews suggest that soil nitrogen (N) dynamics plays a key role in driving 
grazer effects on soil carbon dynamics (Derner & Schuman, 2007; Pineiro et al, 2010). 
Thus, it may be surprising that the model performed well without incorporating any 
explicit N dynamics. One reason for this is that, across the rainfall gradient in the 
Serengeti, rainfall and soil moisture are the primary drivers of productivity, with soil N 
having a much more muted effect (McNaughton, 1985; Holdo et al, 2007). Second, the 
adjustment of productivity for soil texture and water holding capacity in the model 
functions (see Eq. (1)) may have incorporated hidden effects of soil N on carbon 
dynamics, to the extent that soil N is positively correlated with clay and silt (Ruess & 
Seagle, 1994) and grazing intensity (Augustine & McNaughton, 2006). Therefore it seems 
unlikely that the model fit might be improved by an explicit incorporation of soil N as a 
variable governing production and microbial abundance and respiration. 

Another key factor that is important in most other soil carbon models is temperature, 
and it perhaps surprising that the model performed well without a response to 
temperature. In tropical systems like Serengeti, soil temperatures vary from 23° C in the 
wet season when soils are cooled by evaporation to 26° C in the dry season. This low 
temperature range maybe typical for tropical grasslands (Caquet et al, 2012) and is well 
within the normal active range for many soil microbes. The large variation in rainfall 
across sites and accompanying very small variation in mean temperature suggests that the 
principal variable driving soil microbial activity in this system, and perhaps in most 



Ritchie (2014), PeerJ, 1 0.771 7/peerj.233 



] 19/27 



PeerJ 



tropical grassland systems, maybe soil moisture (Caquet et ah, 2012; Manzoni, Schimel & 
Porporato, 2012). It is possible that grazing modified soil temperatures more than was 
assumed here, and feedbacks between grazing, temperature, soil moisture evaporation, 
and microbial activity might improve the predictability of the model. However, including 
such dynamics would appear unlikely, from the good fit of the model to data in this paper 
(Fig. 6) and from previous results ( Wilsey et ah, 2002; Epron et ah, 2004; Caquet et ah, 
2012), to provide a significant improvement in model predictive power for tropical 
grasslands. 

The mass-balance (SNAP) model developed here contrasts in several ways from 
current standard models of grazing impacts on SOC, such as CENTURY (Paustian, 
Parton & Persson, 1992; Schimel et ah, 1994; Parton et ah, 1995) and the Hurley Pasture 
model (Arah et ah, 1997; Thornley & Cannell, 1997). First, these previous models are 
much more complex, as they require many ( > 20) parameters, such as pH, exchangeable 
bases, temperature, soil nutrients, plant tissue nutrients, etc., beyond the five parameters 
required by the SNAP model. The more complex models also consider changes in 
parameter values over relatively fine time scales (days, weeks, or months), as opposed to 
the annual pulse-decay cycle of the SNAP model. Moreover, the complex models track 
two or more SOC pools that differ in their turnover of carbon, while SNAP tracks mainly 
changes in pools of lignin and cellulose. While CENTURY and the Hurley Pasture models 
generally predict well changes in temperate grassland SOC related to production and 
decomposition (McCulley et ah, 2005; Pineiro et ah, 2010), they have not been adequately 
tested for tropical systems and likely contain default functional relationships that do not 
apply to the tropics. Moreover, their structure clearly does not include plant 
compensation to grazing, and they are therefore more likely to predict negative impacts of 
grazing on SOC. 

The simplified model and tests presented here suggests that tracking just the stable 
carbon pool, in this case lignin and cellulose, successfully predicted variation in SOC 
stocks across the Serengeti landscape (Fig. 6). This is supported by the recent observation 
(Frank, Pontes & McFarlane, 2012) that the vast majority of SOC in grasslands is highly 
recalcitrant to decomposition, with very slow turnover rates and residence times of 
100-1000 years. Lignin and cellulose, the primary stable C forms produced by plants, 
likely comprise the major input of stable carbon to the soil following decomposition of 
litter or dung (Andrioli & Distel, 2008; Kirschbaum et ah, 2008; Semmartin, Garibaldi & 
Chaneton, 2008). Dynamics of carbon over short time scales and involving fluxes of labile 
carbon forms (Wilsey et ah, 2002) therefore maybe less important in determining the 
long-term dynamics of carbon that dictate carbon stocks. Consequently, the much simpler 
modeling approach used here might be sufficient to explain variation in carbon stocks. 

Work derived from the Serengeti, which is a unique ecosystem with the majority of 
herbivore biomass in migratory large mammal grazers (wildebeest and zebra), is 
sometimes viewed as inapplicable to other more common grazing systems, or even other 
ecosystems in Africa. While the Serengeti certainly does not exactly mimic most 
livestock-dominated grazing systems, it would be unfair to discount the applicability of 
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Serengeti data and functional relationships to other ecosystems. First, the Serengeti and 
its history of grazing research offer the only such data for tropical grasslands for some 
measures (e.g., root biomass versus rainfall, soil microbial respiration rates versus soil 
properties, etc.). Secondly, from Table 4, it is apparent that the model was tested across a 
range of soil types (sandy to silty clay), mean annual rainfall (450-900 mm), and intensity 
of grazing by both resident non-migratory and migratory grazers. In fact, the site with the 
highest mean grazing intensity features abundant resident grazers that consistently 
remove biomass and maintain a grazing lawn (Anderson et al, 2010), similar to many 
livestock systems. Nevertheless, additional tests clearly are needed from resident wildlife 
grazer systems and free-ranging livestock systems to test the broader applicability of the 
model and specifically the importance of plant compensation in soil carbon 
dynamics. 

CONCLUSION 

Compensation to grazing by plants in Serengeti National Park occurred across a wide 
range of grazing intensities, maintaining LAI at the expense of stem biomass up to 
grazing intensities near 70%. Incorporation of this LAI response into a simplified model 
of soil carbon dynamics, developed with data from the extensive research history in 
Serengeti, suggests that plant compensation can allow grazers to increase SOC stocks at 
intermediate grazing intensities and that this influence is stronger at higher fire 
frequencies. The model performed very well in predicting current SOC stocks at both a 
plot and site level. The success of the model indicates the potential utility of tracking the 
fate of stable carbon (e.g., lignin and cellulose) rather than focusing on fluxes in 
understanding SOC dynamics. 

The predictions of the model, which would likely apply mainly to tropical grasslands, 
suggest that two factors commonly under management control, grazing and fire, can have 
large impacts on soil carbon stocks. In particular, there maybe potential for land users to 
(1) reduce fire in lightly grazed, fire-prone areas of Africa and Australia (Kirschbaum 
et al, 2008; Williams et al, 2008; Ryan, Williams & Grace, 2011; Woollen, Ryan & 
Williams, 2012), and (2) to restore degraded soils in overgrazed areas by better 
management of livestock rather than de-stocking (Conant et al, 2001; Lai et al, 2007; 
Schipper et al, 2007). Such practices could help rebuild organic matter and soil carbon 
across vast areas of the developing world, as tropical grasslands and savannas account for 
at least 10% of the earths land surface and occur primarily in developing countries. The 
simplified model presented here requires only a few measured parameters, so predicting 
changes in accordance with altered land use may be more feasible than using currently 
popular models, like CENTURY, with large numbers of parameters. The model presented 
here (Fig. 1), with its inclusion of plant compensation to grazing, may help explain why 
grazing in C4 grass-dominated grasslands can help sequester soil carbon (Derner & 
Schuman, 2007; McSherry & Ritchie, 2013) and may advance our understanding of soil 
carbon dynamics in tropical grasslands. 
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